Field Theoretical Approach to Electrochemical Deposition* 
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In this work we present an application of the A</>* field theoretical model to the adsorption of 
atoms and molecules on metallic surfaces - the electrochemical deposition. The usual approach to 
this system consists in the computational simulation using Monte Carlo techniques of an effective 
lattice-gas Hamiltonian. We construct an effective model towards a comparison between the lattice- 
^ gas Hamiltonian and the discrete version of the Ai^* Hamiltonian, obtaining the relationships between 

. the model parameters and electrochemical quantities. The model is studied in the mean field 

' approximation, and the results are fitted and compared to numerical simulated and experimental 

data. 

o 

^ , I. INTRODUCTION 

There is a recent convergence on the subjects of studies and techniques of electrochemistry and surface science. 
In particular, the processes where ions and molecules are deposited over a metallic surface can be treated with both 
analytical or numerical statistical mechanics techniques. These processes are of fundamental interest since it presents 
detailed aspects of the electrode-electrolyte interface structure. Electrical current and potential for electrochemical 
deposition in metallic surfaces are obtained by cyclic voltammetry, showing current peaks in a determined potential 
range. The shape, position and number of these peaks depends on the surface where the deposition occurs, as well as 
^ on the electrolyte nature. The peaks are related to phase transitions of the adsorbed layers ( [1], [2]), and this aspect 
O makes this electrochemical subject be well adapted to statistical mechanics studies. 

I ^ I ' The dynamics of the adsorption of atoms and molecules on metallic surfaces in electrolytic solutions due to voltage 
' differences - the electrochemical deposition - can be described, in an effective way, by the lattice-gas model. In this 
model, the adsorption of an element (adsorbed specie) changes the energy of the system due the bound energy of 
the element to the surface, and due lateral interactions between adsorbed species. The system is described by a 
Hamiltonian that reproduces interaction between nearest neighbors (but not only first nearest neighbors) and terms 
\ relative to an external field (adsorbed specie electrochemical's potential). The Hamiltonian is usually treated by using 
■ computational simulation, such as Monte Carlo simulations, in both equilibrium ( [3]- [6]) and quasi-equilibrium situ- 
' ations( [7]- [9]), reproducing surface coverage and voltammetric current profiles in good agreement with experimental 
data. Also analytical treatment is employed in these systems ( [10]- [16]). However, in contrast with the first ones, 
these treatments usually employ some empiric approximations, 
"j^ ■ The analytical treatment of this model is quite difficult, since it involves up to fifth neighbors interactions [9], 
' characterizing the non-local nature of the model. Another aspect that increases the complexity of these treatments 
I ' arises when we consider potential terms of order higher than quadratic. In fact, those higher order potential terms 
' O ' are usually avoided. 

^ ■ The continuous limit of the discrete lattice-gas Hamiltonian is a field theoretical model. In the treatment of such 
model, there's no great difficulties, although requiring perturbative expansions, in dealing with higher order potential 
terms. Nevertheless, non-local terms in field theoretical models are of difficult treatment, and as long as it can be 
done, one looks always for local field theoretical models. 

In this work we present a field theoretical model constructed from the continuous version of the electrochemical 



H \ deposition lattice-gas Hamiltonian. To avoid non-local terms, we will consider only first nearest neighbors interactions, 
. - 1 but, in order to reproduce the phase transition presented by the electrochemical deposition, we will introduce a fourth 
order potential term, leading to the A(/)^ model. 

We must remark, nevertheless, as indicated in [17], that the numerical Monte Carlo treatment of such systems 
has some advantages in the detailed description of the surfaces, of the distribution and sizes of deposited ions and 
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molecules, and so on. But an analytical, or specifically, a field theoretical description of the electrochemical deposition 
process can improve the understanding of the phenomena involved, since it provides simple analytical expressions for 
some of the observable quantities and gives a comprehensible description of how the symmetry breaking is related to 
the phase transition of the system. It also opens the possibilities for the study of the presence of topological solutions 
(solitons and vortices) in the distribution of the adsorbed specie over the electrode surface, which, in future, could 
represent the possibility of new technological features. 



II. CONTINUOUS VERSION OF THE LATTICE-GAS MODEL 



Let us start from the grand-canonical lattice-gas Hamiltonian ( [4], [5], [9]) for the electrochemical deposition of 
one specie in a square lattice, given by 

N N 

H = ~J ^ {cijCij+i + CjjCj+ij} + M ^ Cij + H3, (1) 

where J is the lateral interaction between species occupying nearest neighbors sites, Cij = {0, 1} denotes the occupation 
of sites labeled by its row and column numbers i and j, fi is the electrochemical potential and H3 is higher order 
interaction terms. The electrochemical potential fi is related to the electrode potential through: 

C 

M = Mo + A:Tln— -67!/, (2) 

where Co is some reference value for the specie concentration in solution, /ig is the electrochemical potential at 
the reference concentration and zero electrode potential V, C is the specie concentration in solution and 7 is the 
electrovalence of the specie. 

The thermodynamic conjugate to the electrochemical potential is the surface coverage O, that represents the fraction 
of occupied sites 

N 



e = 7v-i^ci (3) 



i=0 

, where is the total number of unity cells in the electrode surface. The charge transported through the sub- 
strate/electrode interface in adsorption/desorption process is given by 

q = -ejS, (4) 

and the voltammetric current density (voltammetric current per adsorption site) can be obtained from its derivative 

. _dq __ de __ de djidV 

^'"di~~^'^~dt ~~^^diidV~dt' ^' 

where we express cq. (5) as a function of the; scan rate 

Rewriting (1) by using a symmetric occupation number given by 

(f>hj = ^(^iJ - 1; = {-!.+!} (6) 

we find 

N N 



J 

where we had made use of periodic boundary conditions. Introducing the interaction strength by unitary cell 
where o is the lattice spacing and a? is the unitary cell area, we find 

J N /2 ^ 

H = --^Y1 ^^{'t>i,j^i,j+i + ^i,j(f>i+i,j} H ^^^i^^ + ^3 (8) 

«J=1 »)J=1 
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or 



J ^ 1 11 1 



where z is the number of nearest neighbors to the site 

Applying the periodic boundary conditions, eq. (9) reduces to 




(10) 

u / I 

-g^E'^'C-^^E^X.+^a. 

In order to obtain a field model that corresponds to (10), we must assume that the lattice spacing is negligible when 
compared to the lattice dimensions, and that the occupation number of neighbors sites changes slowly from one site 
to another. In fact, as stated by eq. (6), (^^.j can assume only one of two discrete values, ±1. To obtain a genuine 
field theory, we must replace the occupation number ^^.j in eq. (10) by its configurational moan (p{x,y) <pij >, 
where the brackets represents the mean over all configurations allowed to the system at certain fixed conditions 
(temperature, solute concentration, electrode potential). Since cyclic voltammctry is carried out with slow scan rates, 
in a quasi-statical procedure, it seems to be reasonable that all allowed states with some fixed energy could contribute 
to a specific current measurement. 

Denoting H and (p for the configurational mean hamiltonian and occupation number respectively from now on, we 
can take the limit a ^ 0, and replace the lattice columns and rows labels i and j by the sites positions x and y, 
resulting in 

H = l dxdy^^^ (^</.)' - - ~'^+//% | +Hz = l dxdy h, (11) 

where (f) = (f){x,y). Identifying the Hamiltonian density h, we found 

where /13 is the Hamiltonian density of the higher order interaction terms, and 

H3= [dxdyhs. (13) 



By applying the appropriate Legendre transformation, we obtain the Lagrangian density for the continuous model 



Here, we are also neglecting higher order correlations, i.e., we are assuming < <j)^ >=< (j) 
Rewriting the Lagrangian density in terms of 



^=^^ (15) 



we obtain 



Finally, defining the quantities m and c by 



(16) 
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m2 = -^, (17) 

and 



2J-U 

we obtain, for the Lagrangian density of the continuous extension of the lattice gas Hamiltonian applied to the 
electrochemical deposition, 

« = -^('^^)'-^mV+c(/'-/i3, (19) 

formally identical to a scalar field model with a symmetry breaking term. The quantities appearing in (19) are related 
to the lattice gas model effective quantities through eq. (6), (15), (17) and (18). 

The surface coverage, given by eq. (3), can be given in this continuous version of the lattice-gas model by 

" 1 /• 1 

e = iV-i ^ = TV- V- ^ + - = iV-ia- J dx dy cp{x, y) + -, (20) 

where we have made use of (6) and (15). 

The model described by (19) corresponds to the non-linear sigma model, since it is restricted, by means of eq. (6), 
to —1 << <p >< 1, wich implies in 

In this work we are interested in the description of the phase transition in electrochemical deposition, which occurs 
in the region far from saturation (tp = ±-\/j /2). Thus, we will treat the as a non-restricted field, keeping in mind 
that the results will be valid only in the region where |<^| < VJ/2. 



III. THE \ip* MODEL IN THE MEAN FIELD APPROXIMATION 

From eq. (19) we can see that, with vanishing breaking term and higher order interactions, the Lagrangian for the 
continuous version of the lattice gas model presents a symmetry related to the exchange of the occupation number 
(fi —cp. Since the transition between the state where the adsorption of the specie is favorated and the state where 
desorption is favoured is dictated by the electrode potential, present in the explicit symmetry breaking term (18), 
we assume that the higher order interaction terms will be also symmetric to the exchange of the occupation number, 
and thus only even order terms must be present in /13. The lowest even order interaction term is the (fi^, and we will 
assume 

hs = (22) 

where A is the coupling constant and was inserted here for dimensionality reasons. By this relationship, we are 
assuming that the potential energy of the system is non-quadratic in ip. Prom (19) and (22) we obtain 



I = 



|(^^)'-^mV' + c(^-^^', (23) 



the A(p^ model with an explicit symmetry breaking term. Due the fact that the mass term in (23) is positive, since 
is negative (17), the model presents spontaneous symmetry break in the limit c = 0. In this limit, the minimum 
of the potential energy in (23) occurs in 



Vn,in = \ — T =V (24) 



or 



4 



A 



Perturbativc solutions must be constructed around one of these minima. Introducing the shifted field = ■ 
in the Lagrangiau (23), we find 



2 '>t'\\ 77A 3 A 



"^^"3!^; - 3! '^^ - 4!^^- 

Lagrangian (25) is no longer symmetric in the sense that it is not invariant under the exchange ips = —fs, so the 
model shows an spontaneous symmetry breakdown. The potential energy minimum occurs now at the point = 0. 
Applying the EuUer-Lagrange equations 

^( dl \ dl 

= (26) 



,2 



in (25), we obtain, for the dynamic equation of the system. 

We can proceed the solution of eq. (27) perturbativelly around tps = 0. Thus, expanding ips as a power series of a 
pcrturbative parameter e, i.e., = Vo + £<Pi + e^V2 + with (po = 0, and keeping only the lowest order terms {e^) 
implies in 



2 

— Tn 

6a^ 



m^r] + c = (28) 



, a cubic equation that can be solved analytically for 77. As showed in figure 1, equation (28) presents one (figures la 
I and III) or three real solutions (figure la II), for the potential energy as a function of (ps, corresponding to the local 
minimum or maximum. Stability is granted for the global minimum, for c = we have two symmetric global minima, 
but for c ^ the system has this symmetry explicitly broken, as one can see from figures lb - Id. Thus, in c = the 
system presents a transition between one global minimum to another. 
Solutions to equation (28) is given analj^ically by 



and 



with 



^3=(-^ + ^ -V3.^ + ^), (31) 



A = 3ca^A2 + AV ^ . (32) 

We must observe that solutions (29), (30) and (31) are independent of the unitary cell area, since it appears only 
in the combinations m^a^ and ca^, and, due (17) and (18) 
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m^a? = — = -z (33) 



and 



2 2J — jJL 2 2J — fj, 



Thus, the number of free parameters of the model is reduced for only four: the electrochemical potential /x, the 

effective lateral interaction J, the strength of the non-quadratic local interaction A and the effective electrovalcnce of 
the adsorbed specie 7. The latter is also limited into a range of values near the electrovalcnce of the adsorbed specie. 
First order (s^) solutions to (27) implies in 

VVi = (^+m2)vi. (35) 

Solutions of eq. (35) can improve the results presented in this work, with oscilatory solutions around the mean 
solution ipo. 

IV. SOFT BREAK LIMIT 

In the limit of vanishing c, solutions to eq. (28) reduces to 



Vc=o = 0; r/c=o = ±\l ^ • (36) 

For very small c, we must suppose ij = r]c=o + S, where 5 is a small perturbation around the c = solution. Replacing 
these perturbative solution in eq.(28) and retaining only first order terms in S, one can find 



^- = ±V^^ — ^- ^''^ 

This approximate solution for t] allows us to obtain some simple estimative expressions for the voltammetric current 
and STirfacc coverage that guide ns in the initial parameters choice for the non-linear fitting of the model to the 
experimental data. The physical measurable parameters are the electrode potential V, the scan rate dV/dt, the 
voltammetric current i and the surface coverage 0. 

The energy of the system in the state represented by solution (37) can be read from eq. (25). The solution that 
represents the global minimum energy switches from rj^ to rj- when c runs from negative to positive values. The 
electrode potential at the transition point, obtained from (18) and (2) for c = is given by 

K = Z1£±J^, (38) 

67 

where we call /Uqc = Mo + kThi{C/Co), the electrochemical potential at solute concentration C and temperature T. 
The current at the transition point can be obtained from (5), where the intermediate derivatives can be obtained from 
(18) and (2), and, from (28), we have 

— ^ri'^dr] - m^dri + dc = (39) 
2a^ 



I.e. 

dr} 



dc Xm'^la^ + 2m2 ' 

Replacing these results in (5), we obtain 



(40) 



i = £V I ^ (41) 
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a result that is not restricted to the soft break hmit. 

At the transition point, where rj is given by the non zero solutions of (36), replacing (17), we obtain 



Near c = 0, with 77 given by (36), is quadractic in V, as one can see from (41), (37), (18) and (2). The critical 
electrode potential can be evaluated by means of eq. (2) and (18) with c = 0. 

The fact that is quadractic in V near the critical point, together with cq. (38) and (42) can be used to set the 
initial guess parameters for a non-linear fitting of eq. (41) to experimental data. 



V. NUMERICAL RESULTS 



Voltammctric current and surface coverage for the process of clcctrodcposition of one single adsorbed specie on a 
square-lattice metallic surface at constant temperature and solute concentration can be reproduced in lowest order 
of the mean field approximation by equations (41) and (20) . The free parameters of the model are J, /Uoc, A and 
7. The latter one is restricted to values near the electrovalence of the adsorbed specie. In comparison with the usual 
Monte Carlo approach, it is a reduced set of parameters, since in this approach, besides the parameters 7 and Hoc, 
up to five lateral interaction between species is used in general [9]. 

Figure 2 shows the results of the fitting of eq. (41), with 77 given by (29) to (31), to experimental and simulated 
data for the electrosorption of Br in the surface of A(/(100) electrodes ( [18], [19], [20]). We use a non-hnear least 
squares fitting, applying the steepest descent procedure. Initial values for the parameters are estimated using eq. (38) 
and (42), and the Br electrovalence. The best fit results in A = —47.7 x 10'', /igc = 610meV and J = AA.QmeV. 
AUthough J and /xqc are strongly dependent on 7, the set of parameters gives a fit of the ratio i/{dV/dt) that is 
independent of 7. The surface coverage, neverthelles, depends on the value of the effective electovalence, and it was 
used to fix 7 = —0.73, which corresponds to the value given in [18] and gives a good agreement between the model 
and experimental data [19]. This value is also near the effective electrovalence estimated as 7 = —0.71 in [21], and no 
significative differences are found in the results computed with 7 = —0.73 or 7 = —0.71. Experimental CVs present a 
shoulder followed by a sharp peak, this second corresponds to a phase transition from a disordered to an ordered phase 
[19]. The shoulder is atributed to configurational flutuations on the Br adsorbed layer [18], and of course cannot be 
reproduced in the model. Due the size of adsorbed bromine being bigger than the ^4^(100) lattice spacing, in the 
ordered phase the bromine is distributed in an arrangement that covers only 50% of the Ag surface. Results for the 
surface coverage presented in figure 3 arc normalized to reproduce this fact, and shows an reasonable agreement with 
the experimental data. 

Figure 4 shows the results of the fitting of the model to experimental data for the Urea over Pt(lOO) surface [22]. 
Fitted parameters are A = —58.4 x 10^, Hoc = 97.8meV, 7 = — 1 and J — 3.04meV" . Figure 5 shows the surface 
coverage for this process. In figure 4 the dotted line corresponds to the quadratic fitting of We can see that there 
is an exceptional agreement between experimental data and the quadratic fitting in the sub critical potcmtial region. 
Of course, this result is beyond the model presented here, since this quadratic solution do no take into account higher 
order terms in c on the solution of eq. (28). 

On figure 5 the dashed horizontal lines correspond to the limits < tp >= rj = /2, values of 77 inside this region 
and relatively far from these limits corresponds to the linear approximation of (19). As one can see from figures 
(4) and (5), the phase transition region is inside these limits. In particular, for the Br adsorption in Ag, all the 
voltammogram depicted on figure 2 is inside the region that corresponds to the linear approach. 

Both results shows a good agreement between calculated and experimental/simulated data. The computational 
advantage of the field theoretical approach over the Monte Carlo simulation is the reduced number of parameters of 
the former. 



VI. CONCLUSION 



In conclusion, we constructed a field theoretical model to describe electrochemical deposition. Based on the lattice- 
gas hamiltonian of the process, and taking a configurational mean of the occupation numbers, we obtained, in the 
continuous limit, a scalar field hamiltonian whose coefficients are related to the lattice-gas model effective quantities. 
Including higher order potential and treating the model in the mean field approximation, we obtain a single system of 
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equations ((20), (29)-(31) and (41)) that reproduces the current versus electrode potential behaviour in electrochemical 
deposition. The fitting of the model with experimental and numeric simulated data was obtained and shows a good 
agreement between model and experimental/simulated results. In comparison with the Monte Carlo approach, the 
model presented here has the advantage of presenting a reduced set of free parameters and simple analytical expressions 
for the observed quantities. 

These results open the possibilities of new analytical description of other more complex deposition process, such as 

two species deposition, where one must use a two scalar field model, introducing new symmetries in the model (work 
in progress). Other even order potential terms could also be used to reproduce phase transitions between different 
ordered phases. They also suggests that the non-linear description of the process can presents other solutions, such 
as topological solutions like solitons and vortices. In fact, the Xcp'^ term can be viewed as a fourth order term in the 
expansion of a cos((/)) term in the lagrangian (19), leading to the sine-gordon model [23], which presents solitons. 

Finally, the hamiltonian studied in this work shows us that the phase transition presented by electrochemical depo- 
sition is related to an explicitly symmetry breaking process, dictated by an external field (electrochemical potential), 
controlled by the electrode potential. 



VII. ACKNOWLEDGMENTS 



The authors are grateful to CAPES-Brazil. We also thanks A. R. Pereira and M. C. Nemes for their suggestions to 
the present work. 



[1] J. G. Dash and J. Ruvalds, Phase Transition in Surface Physics, (Plenum, New York, 1980). 

[2] D. Nicholson and N. G. Parsonage, Computer Simulation and the Statistical Mechanics of Adsorption, (Academic Press, 

New York, 1982). 
[3] D.P. Landau, Phys. Rev. B 27, 5604 (1983). 

[4] P.A. Rikvold, J.B. Collins, G.D. Hansen and J.D. Gunton, Surf. Sci. 203, 500 (1988). 
[5] J.B. Collins, P. Sacramento, P.A. Rikvold and J.D. Gunton, Surf. Sci. 221, 277 (1989). 
[6] P.A. Rikvold, M.R. Deakin, Surf. Sci. 294, 180 (1991). 
[7] M. Gamboa-Aldeco ct. al.. Surf. Sci. Lett. 221, 297 (1993). 

[8] P.A. Rikvold, M. Gamboa-Aldeco, J. Zhang, M. Han et. al.. Surf. Sci. 335, 389 (1995). 

[9] P.A. Rikvold, J. Zhang, Y.E. Sung, A. Wieckowiski, Coll. Surf. A 134, 3 (1998). 
[10] D.A. Huckaby and L. Blum, J. Chem. Phys. 92, 2646 (1990). 
[11] L. Blum and D.A. Huckaby, J. Chem. Phys. 94, 6887 (1991). 
[12] D.A. Huckaby and L. Blum, J. Electroanal. Chem. 315, 255 (1991). 
[13] D.A. Huckaby and L. Blum, J. Electrochem. Sac. Conf. Proc. Ser.92-1, 139 (1992). 
[14] L. Blum and D.A. Huckaby, J. Electrochem. Soc. Conf Proc. Ser.93-5, 139 (1993). 
[15] L. Blum and D.A. Huckaby, J. Electroanal. Chem.375, 69 (1994). 
[16] L. Blum, D.A. Huckaby and M. Legaut, Electrochem. Acta 41, 2207 (1996). 
[17] G. Brown, P.A. Rikvold, S.J. Mitchell and M.A. Novotny, cond-mat/9805126 
[18] S. Mitchell, G. Brown and P. A. Rikvold, Surf Sci. 471, 125 (2001). 
[19] B.M. Ocko, J.X. Wang and T. Wandlowski, Phys Rev Lett79,1511 (1997). 
[20] G. Valctte, A. Hamelin and R. Parsons, Z Phys Chem 113, 71 (1978). 

[21] I. Abou Hamad, Th. Wandlowski, G. Brown and P. A. Rikvold, J. Electroanal. Chem. 554-555, 211 (2003). 
[22] P.A. Rikvold, J. Zhang, Y.-E. Sung and A. Wieckowski, Electrochem. Acta 14, 2184 (1996).. 
[23] G.B. Whitham, Linear and Nonlinear Waves, Wiley, 1974. 



8 



Figure Captions 



Figure 1: (a) Cubic equation corresponding for the configurational mean value of the occupation number (f). The 
horizontal dashed lines correspond to the 77 = axis with I) c< 0, II) c=0 and III) c> 0. Note that for c> and c< 
one can have one or three real solutions, (b), (c) and (d) Potential energy as a function of rj for different values of 
the symmetry breaking term. For c=0 the two minima are symmetric, as c goes from negative to positive values, the 
global minimum switches from one side to another, corresponding to the phase transition of the system. 

Figure2: Current density as a function of the electrode potential for the electrodeposition of Br on ^^(100) in the 
(j)^ model (solid line). Full and empty squares correspond to simulated results [18], empty squares are points not took 
into account in the nonlinear fitting. 

Figure 3: 0^ model results (solid line) and experimental results (full squares) [19] for the surface coverage in the 
electrodeposition of Br on Ag{lQQi). 

Figure 4: Current density as a function of the electrode potential for the electrodeposition of Urea on P^(IOO) in 
the (/)^ model (solid lino). Experimental results (full squares) are obtained from [22]. Dotted line corresponds to a 
quadratic polynomial fitting oH~^ . 

Figure 5: Surface coverage for the the electrodeposition of Urea on Pf(lOO). Dashed horizontal lines correspond to 
the approximate limits of the linear approach presented in this work. 
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